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To check the reliability of the Majorana representation of quantum spin systems, we use it to 
compute the ground state energy of an antiferromagnetic spin 1/2 chain to one loop order. We find 
a very small one loop correction of the mean field energy and a discrepancy of > 100% compared 
to the Bethe Ansatz result. We conclude that a careful handling of the gauge degrees of freedom of 
this representation is crucial to get correct results. 
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I. INTRODUCTION 

IS ' 

( \ The Majorana representation is a promising tool for describing frustrated quantum spin 1/2 antiferromagnets and 
heavy fermion systems jl[ , . It provides a description of spin operators in terms of bilinears of real fermions without 
any constraints. It is useful to consider the representation of spin operators of the group O(N) and take N = 3 only 
C$ . at the end: 

O(N) :s k J = ~ (r£r£ - rjlrt) , k, I = 1, 2, ..N, [rfc, rf y ) + = S kl S xy , rj k x + = r, k x 
§ ! 0(3) : sZ = -~r%XriZ (1) 

A mean field treatment of spin systems in this representation gives rise to O(N) invariant bond variables in a natural 
way and, therefore, this representation seems suited for a description of the disordered or "liquid" phases of quantum 
spin systems. Indeed it is not difficult to see, by using the Majorana representation, that most O(N) spin systems 
(in the lowest representation that generalizes spin i) reduce to RVB like spin liquids at N = 00 Q]. However, for the 
Majorana representation to be really useful, it must also be able to describe the competition between ordered and 
I liquid spin systems as a function of the amount of frustration and for N = 3. 
7-H ■ Recently, it was pointed out that a simple Hartree Fock treatment of the spin 1/2 Heisenberg antiferromagnetic 
chain gives qualitatively correct results [0. The purpose of the present note is to show that in a more systematic 
approach the one loop correction to the mean field energy is small and that the ground state energy still differs from 
the exact Bethe Ansatz result by more than 100%. 

It is very likely that the local Zi invariance of the Majorana representation is responsible for this discrepancy and 
that it must be dealt with correctly in order to obtain reliable results. The situation is reminiscent of the Yang Mills 
theory of strong interactions in which a correct treatment of the local invariance is crucial. 

Although the local Zi invariance of the Majorana representation was previously dealt with in the context of the two 
channel Kondo problem || a more systematic treatment of this invariance using standard gauge theory techniques is 
O ■ needed. 
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II. MEAN FIELD 



The Hamiltonian of the O(N) antiferromagnet is 



o(n) .h= j2 r k J4 l -j = \ E (») 2 

<x,y>k,l <x,y> 

0(3) :H= £ -* a . = ! E (t,-^) 2 (2) 

<x,y> <x,y> 

with < x, y > denoting neighboring points, J = 1 and we will suppress the internal symmetry indices in the following. 
To allow for a systematic treatment of fluctuations the mean field must be introduced by a Hubbard-Stratonovich 
transformation: 
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Z = I re- 1311 = J Drje'fo dt (i^v^+H(r,)) 

j DBDne-So^ 9 ^ 11 ^, H(B, rf) = W^ + \Y, B ^ 

J ^ J xy ^ 

x,y u x.y 

Z = [ DBe~P F ^ , F(B) = - -L log det (d T + iB) (3) 
J 2p 



For a d = 1 chain and with a uniform auxiliary field 



M 

B X y = (Sx,y + 1 ~ Sy„x+l) (4) 



we find 



H(B, 17) = -M 2 - yj] (VphV-pi* - iI-phVph) sin? 

= <W,o (5) 

where L denotes the length of the chain. The mean field saddle point is M = — and gives a quasiparticle dispersion 
and ground state energy (for the original Heisenberg Hamiltonian) of 

E(p) — sinp ^= 3 — sinp (Bethe Ansatz: — sinp) 

7T 7T 2 

< H ^> + E = + H ™ _ .08 (Bethe Ansatz: - 0.44) (6) 

L 8 2-k z 8 

It is known that the optimal mean field in the Majorana representation is completely dimerized for a large class 
of Heisenberg antiferromagnets in any dimension (J] and we could check whether the fluctuations about the mean 
field restore the uniform mean field configuration in d = 1. However, for our purpose of demonstrating a bug in an 
uncritical use of the Majorana representation this is unnecessary since: 

• the exact Bethe Ansatz solution does not break translations - 

• as we will see the fluctuations about the uniform mean field in the Majorana representation give only a small 
correction to the mean field energy, with the one loop corrected ground state energy still differing from the exact 
result by more than 100%. 

III. ONE LOOP CORRECTION 

The self consistent Gaussian approximation is a very reasonable approximation that gives the first 3 terms of the 
Stirling series of T(n) = J* °° dxx n ~ 1 e~ x at one stroke and which should, a priori, give good results also for the spin 
chain in the Majorana representation. In this approximation, the partition function reduces to 



e-™ _ 1 . . . S 2 F_ 

2 



e-W, F(B) = F(B) + - log det (7) 

det 2 P SB 



TW 1 

To identify the quadratic form |^one must expand the free energy F(B) about the saddle point 

F(B + SB) = F{B) + l^SBSB + O (SB, SB 3 ) 
2 SB Z ' 

klT5S 5B5B = T3 f ' dT Y, SB ly+T3 tr (jT^5) SB (V^B ) SB 
2 SB 2 4(3 J ^ v A/3 \d t +iBj \d t +iB J 

After some simplifications the total quadratic form turns out to be 



(8) 
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6 ^6B6B~^6b*(q)F(q)8b(q) 

F ^ = 1 + TS £ (l-e i(r - p) )ff(pMr),?=(9o,gi) 

p-4-r — q— 
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5B = -^^«(g)e^, 

^ V g{p)e ipx , g(p) = — (9) 



dt+iB y/Lj3^ ipo + Msinpi 

At zero temperature, the sum over Matsubara frequencies turns into an integral and one obtains 

1 fS 2 F\ T =a 1 f°° r 

5E f i uctuations = — log det I \ = -^J dq J dq 1 \ogF(q () ,q 1 ) 

F(qo, qi) can be expressed in terms of elementary functions: 

l v , 2lp , n F (Msm(p- f))-n F (Msin(p+f )) 
*W> qi) L^ 1 G > Msin(p-f)-Msin(p+f)-i go 

— 2 

1 / 2 sin pcosp 



(10) 



dp 



(11) 



1 / Vl + x 2 ^ ! / singi/2 \ \ go 

= 1 + ^^ 1 + ^(fy tanh Ittt^ J J ' x = ^hTf 

The total energy that includes the fluctuations about the un dimcrizcd saddle point is, therefore, 

_ M 2 NM 

En — ~Z n (" OE fluctuations 

8 2tt 
M / ,7r r°° q 
SEf luctuaUons (M) = I dq dx log F sin -dxdq (12) 
Jo Jo 2 

Setting N = 3 and M = — we find SEfi uctua tions = 0.036 and since the correction is so small, it is not necessary to 
do the Gaussian approximation self consistently. The total energy differs from the exact Bethe Ansatz result still by 
much more than 100%. So the loop expansion converges very well, but. ..to the wrong result. 

IV. THE PROBLEM OF GAUGE INVARIANCE 

In the Majorana representation of eq(l) , there is an infinity of operators that commute with the Hamiltonian pi 



N 



Q(li) = i]] »?sp 

u=l 

[Q( Vx ),H} = 

N even : [Q(r? x ), Q y (r) v )]_ = 

N odd : [Q( Vx ),Q y (r) v )] + =^~\ 5 x , y (13) 

The oscillation in statistics of Q(r) x ) as a function of N may be considered to be responsible for the even/odd 
oscillation in the energy of a small cluster as a function of N that was observed in j| . The r\ x are gauge transformed 
by commuting (anticommuting) them with the operators Q{"q x ). The theory of quantising gauge invariant systems is 
highly developed but a caricature of existing methods runs as follows: 

• some redundant parameter is eliminated and the system, now without gauge invariance, is treated in terms of 
the remaining variables by conventional methods - 
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• or, alternatively, one refrains from eliminating redundant variables. Instead, one imposes a constraint on the 
space of variables and provides a correction factor in a way that essentially extracts the volume of the gauge 
group from the partition function. This is done by inserting unity in the form of 

1 = A FP (V) I De6(f( Ve )) (14) 



into the partition function, where f(rj) is a function of r\ that breaks the symmetry, and where De integrates 
over the group [gj. App(rj) is the determinant of the operator ^ and is represented by auxiliary or "ghost " 
particles. 

To apply the first method, one notices that the constant multiplicity 2 L / 2 of extra spurious states in a system of L 
spins can be understood in the representation 

Vxfj, — &xfj,fx (15) 

where f x are Majorana fermions and a x ^ are Pauli matrices of the group O(N). It can be shown that the extra 
fermions f x get eliminated by imposing a positive spectrum , irj x rj y > 0, on a non overlapping covering by dimers 
(x,y) or on dimers between points on the original system and an extra copy of it. In the second approach, because 
the symmetry group is fermionic, Afp(i]) is actually the inverse of the determinant of ^ and will be represented by 
bosonic auxiliary or "ghost " particles. The details of this gauge approach to the Majorana representation still remain 
to be worked out, however. 

In conclusion, we have seen that ground state energy of a spin 1/2 antiferromagnetic chain treated by Majorana 
representation and mean field plus Gaussian fluctuations comes out incorrectly, although the loop correction is small. 
This strongly suggests that a proper treatment of the Majorana representation of spin systems must address the 
problem of local Zi invariance or the local supersymmetry of the system. 
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